Magnetic behavior of a non-extensive S-spin system: possible connections to 

manganites 
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We analyzed the magnetic behavior of a S-spin system within framework of the Tsallis nonexten- 
sive statistics, employing the normalized approach. Unusual properties on magnetization, entropy 
and susceptibility emerge, as a consequence of nonextensivity. We further show that the nonex- 
tensive approach can be relevant to the field of manganites, materials which exhibit long-range 
interactions and fractality, two basic ingredients for nonextensivity. Our results are in qualitative 
agreement to experimental data in Lao.67Cao.33Mn03 and Pro.5Cao.5Mno.95Gao.05O3 manganites. 
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I. INTRODUCTION 

In 1988, Tsallis|Q] proposed a g-dependent entropy 
functional which generalized the standard Maxwell- 
Boltzmann definition, to include non-extensive systems. 
In this formalism, "q" is called the entropic index, and 
measures the degree of nonextensivity of the system. A 
few years latter, Curado and Tsallis[Q revised the for- 
malism introducing the unormalized constraint to the 
internal energy and showed that from this entropy func- 
tional a nonextensive thermodynamics could be derived. 
This generalization included the classical thermodynam- 
ics (Maxwell-Boltzmann) , which could be recovered when 
q is set to unity. 

After Curado and Tsallis work, this formalism has been 
successfully applied to various physical systems, where 
Maxwell-Boltzmann framework fails. These include: self- 
gravitating systemsjH |K turbulence H, ^, (t| [^anoma- 
lous diffusion P, |Tc|, [ill p|) velocities of galaxies iful, solar 
neutrinos JU|, etc. 

In spite of these successes, some drawbacks were iden- 
tified in the early formalism, namely: (a) the density op- 
erator was not invariant under a uniform translation of 
energy spectrum; (b) the g-expected value of the iden- 
tity operator was not the unity; (c) energy was not 
conserved In a more recent work, Tsallis et al.]l5| 
circumvent these difficulties introducing the normalized 
constraint to the internal energy of the system. 

Concerning applications to magnetic systems, the un- 
ormalized formalism was used by Portesi et al.Jl(| and 
Nobre et al.jlTj, to describe the paramagnetic behavior 
of a system with N spins 1/2. The authors found a non- 
measurable magnetic susceptibility, since it was exponen- 



tially dependent on the number N of particles. Martinez 
et al.jl8|] analyzed the same S=l/2 system in the frame- 
work of the normalized formalism []l5|], finding a similar 
result to that reported by the Portesi et al. and Nobre 
et al. 

In the present work, we analyzed the paramagnetic 
behavior of N spins S, within the normalized formalism 
and showed that the system effective temperature T does 
not relate with the inverse of the Lagrange multiplier 0, 
as it is normally assumed, but with the inverse of a re- 
scaled parameter (3* (T — l/k/3*), already introduced 
by Tsallis |jl5| . In this approach, the unreal result found 
by Portesi[p"6|| , NobreJlTj] and Martmez[]l8| was no longer 
found and the susceptibility became proportional to num- 
ber of particles, as it is experimentally expected. In these 
early works, no attempt was made to correlate theoretical 
results to experimental magnetic systems of any kind. 

We also suggest that the manganites are physical sys- 
tems where the present concepts can be tested. Experi- 
mental data of Amaral et al.[{L9) and Hebert et al.|Q, are 
in qualitative agreement with the results here reported. 
An analysis of non-extensivity in manganites in the fer- 
romagnetic phase, also taking T — l/k/3* as an effective 
temperature, was already published else where 



II. THE MODEL 

We consider the Hamiltonian of a single spin S in a 
static and homogeneous magnetic field: 



H = -Jx-B = -fi z B = -g[i B S z B 



(1) 



where g=2, for J=S. The g-generalized magnetic mo- 
ment thermal average (ft z ) q is: 
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where (S z ) q is g-generalized thermal average spin oper- 
ator and i?s q) is the generalized Brillouin function. This 
first quantity can be determined in the Tsallis framework 
of normalized ^-expectation values Jlq| , from: 



B. Normalized Tsallis Entropy 

If the entropy functional S q is defined as proposed by 
Rajagopal and AbeKj]: 



Tr{p«S x } 
Tr{pi} 



(3) 



where p is the density operator, derived from the max- 
imization of the entropy. Below, we analyze two different 
proposals for the entropy and show that re-scaling the 
Lagrange parameter, the same density operator emerges 
from both definitions. 



_ Tr{jPS q } _ k 
q Tr{p«} (q - 1) 



1 



Tr{pi} 



(10) 



the density operator p can be determined in a similar 
way as above, yielding: 



1 

p = y 



l-(l-q)PTr{pP}(H-U q ) 



(11) 



A. Tsallis Entropy 

The entropy of the system is defined asjl], 0, [T^ |: 

S q = Tr{p«S q } = Tr{p«}} (4) 



where, 



S q = -khxgp 



(5) 



and k a positive constant. The generalized logarithm 
is defined as 0221: 



In, / = (1 -«?)-*(/ 



-l (f l-q 



1) 



(6) 



The corresponding density operator p can be deter- 
mined from the maximization of the S q functional, sub- 
jected to q-normalized constraint p5||: 



U„ 



Tr{pPH} 
Tr{pi} 



(7) 



and the normalization of density operator Tr{p} = 1, 
leading to: 



!-(!-<?) 



P 



Tr{pi} 



(H - U q ) 



(8) 



where Z q is the generalized partition function: 



Z q = Tr 



l-(l-«) 



P 



Tr{pi} 



(H - U q ) 



(9) 



and is the Lagrange parameter associated to internal 
energy constraint p~5|, In this case, the entropy has a 
well defined concavity, for any value of (?, being concave 
for q > 0, and convex for q < 00 § [fj|, §§. 



Z q = Tr 



l-(l-q)l3Tr{p q }(H-U g ) 



(12) 



In this case, the q parameter is restricted to the interval 
< q < 1, preserving the entropy concavity pi). 

The density operators emerging from both cenarios 
(Eq.S and O) can be re-written in the form: 



1 



Z' q = Tr 



where, 



p* = 



1 - (1 - q)l3*H 
l-(l-q){3*H 

P 



Tr{pi} + (1 - q)U q f3 
for the case A, and: 



0* 



P 



feffj + (1 - q)U 9 p 



(13) 
(14) 

(15) 
(16) 



Tr 



for case B. Here, we suggested that the effective tem- 
perature is: 



T = 



1 

k[3* 



(17) 



and the density operator becomes independent of the 
initial entropy functional. 

The magnetic behavior of a S-spin system will be ana- 
lyzed as a function of the parameter x* = gp, B SB/3* . In 
fact, in terms of a:*, the generalized Brillouin function, is 
given by: 



+S 



B 



(</) 



£ m s [l + (l-q)x*^ 

1 m s = -S 

£ [l + (l-q)x*?f]^ 

77l s — — S 



(18) 
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It is to be remarked that cut-off procedurefl23], |25| , |2q , 
|27j| implies that those states that do not satisfy the con- 
dition: 



l + (l-?Ky>0 



(19) 



must be excluded from the summation. In other words, 
these states are assigned with a zero probability ampli- 
tude, preserving the positive definition of the density 
operator[^6[. 

For a system constituted by N spins S particles the q- 
generalized spin operator thermal average (Eq||) , can be 
written as: 



The occupation probability (OP), as a function of x* , 
for each energy level of S=5/2 are displayed in figure 
3(a) (b)(c), for several values of q. Figure 3(d) shown the 
same quantity for S=2 and q=0.1. We observe that the 
OP does not vanish for negative energy levels (m s > 0), 
even for very large values of x*, in contrast to what hap- 
pens in the case q=l (Maxwell-Boltzmann) . From figure 
3(c) we can see that, for q=0.1, the OP of the positive en- 
ergy states (m s < 0) vanish sharply at the same x* values 
as the kinks observed in figure 1. This occurs for q<0.5. 
Therefore, we correlate the kinks observed on the mag- 
netization curves to the lost of occupation of the most 
energetic states. This is consequence of Tsallis' cut-off. 
In fact, from Eq.|l^, the x* values where the kinks occurs 
can be derived, following: 



(S z )q,N — 



+NS 

E Y(m s )m s [l + (l-q)x*^] 1 - 

m B = -NS 

E Y(m s )[l + (l-q)x*lf]i-« 

m a = -NS 



where, Y(m s ) is the multiplicity and : 



NS 



J2 Y(m s ) = (2S+l) N 



(20) 



(21) 



m s = -NS 



In the particular case of N 1/2-spins particles, Y(m s ) 
is simply given by: 



Y(m s ) 



Nl 



(f-m s )!(f +m s )! 



(22) 



III. RESULTS AND DISCUSSION 

When q is different from unity, general expressions for 
magnetic observables are difficult to calculate and inter- 
pret, even for simple systems. This hinders the compari- 
son between theoretical predictions and experimental re- 
sults. Numerical methods, on the other hand, provide a 
means to directly calculate observables, allowing a deeper 
understating of the theoretical results. 

Figure 1 displays the generalized Brillouin function 
£?s q) vs. x* (Eq.[lq) , for different values of q and S=5/2. 
For q up to 0.5, a series of kinks appear in the curve. One 
also notes that the saturation value -Bs q) |sAT decreases 
with decreasing q. This is more clearly shown in figure 



2(a) and (b), where Bg 



is plotted as a function of q, 



for half- integer and integer spin values, respectively. The 
behavior of a classical spin is also included and can be 
exactly calculated from Eq.[T§with S — > oo, giving: 



kinks 



s 



(1-9) 



(24) 
fl/2 



It is to be remarked that for a half-integer spin, S 
kinks occur, and S for an integer spin. 

Figure 4(a) (b) display, respectively, the unormalized 
(Eq.Q) and normalized (Eq.^) entropy S q , for S=5/2 
and different q values. Note that, for q<0.5, the kinks 
discussed before are present. Besides, for any q<l, the 
entropy does not vanish in the limit x* — > 0. In other 
words, even for high field and/or low temperature, the 
system has a finite entropy, which prevents a fully mag- 
netized state. These features are valid for both entropy 
functionals (section II). 

A general expression for the magnetic susceptiblity can 
be deduced from Eq.|| and Eq.[l^: 



Xq 



lim 



dB 



T 



(25) 



where = qC^ 1 ' is the generalized Curie constant. 
Figure 5 displays Xq 1 as a function of x* for q=1.0 and 
0.8, with S=5/2. 

In order to compare the predictions of the proposed 
model to experimental data, we must investigate how 
the magnetic observables discussed before scales upon 
increasing the number of particles in the system. An 
extensive quantity T scales as 



N 



(26) 



where N is the system number of particle. 

Particularly useful for the purpose of comparison, is 
the magnetic susceptibility. Portesi et al.Jl(| and No- 
bre et al.jf?]], considered the generalized magnetization 
on paramagnetic fase of a N 1/2-spin system in the un- 
ormalized formalism|||, |l|, [jl}, [l?], [23). They found for 

Xq- 



B< q) l 



(23) 



Xq 



T 



(27) 
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where, 



r(q) _ (ffMe) 2 N 



(28) 



and T = l/fc/3 is the temperature of the system. Here, 
P is the usual Lagrange parameter associated to the non- 
normalized constraint of internal energy. 

Therefore, as N increases, the magnetic susceptibility 
becomes infinity for q<l and zero for q>l. In other 
words, there is no linearity among \q an d N. This phe- 
nomenon is called dark magnetism, in analogy to the 
cosmological concept of dark matter [|jl^| and references 
there in]. 

S. Martinez et al.JTq], analyzed the paramagnetic be- 
havior of the same N 1/2-spin system, however, within 
the normalized formalism |l5|, gsj, and found a similar- 
result as Portesi and Nobre for Xq- 

Our proposal is that the paramagnetic behavior of a 
non-extensive N 1/2-spin system should be analyzed in 
the normalized formalism, using the density operator de- 
scribed in Eq.|l3], taking /?*, instead /3, inversely propor- 
tional to the system temperature T. By doing so, the 
generalized paramagnetic susceptibility becomes: 



Xq 



^S = l/2 

T 



(29) 



which does not diverge as N increases, and is indeed 
proportional to N. Besides, this result is independent of 
the choice for the functional entropy. 



IV. POSSIBLE CONNECTION TO 
EXPERIMENTAL RESULTS 

Amaral et al.|l9) discussed the magnetic behavior of 
manganese oxides, namely, Lao.67Cao.33Mn03 and veri- 
fied steps on the curve of M -1 vs. T in the paramagnetic 
phase, as shown in figure 6(a). These steps are analogue 
to those shown on figure 6(b), that represents also the in- 
verse of magnetization as a function of a;* -1 , for q=0.1, 
and encourage the idea of manganites as non-extensive 
objects^]]. The authors relate the change in the slope 
of the curve as an indication of cluster formation, which 
change the effective moment of Mn ions. These clus- 
ters could give rise to fractal structures, as discussed by 
Dagottopg|], and therefore are in accordance to the ideas 
discussed here. 



In addition, Hebert et al.|2Cj found magnetization 
curves in Pro.5Cao.5Mno.95Gao.05O3 that are also in qual- 
itative agreement to those curves in figure 1. 



V. CONCLUSION 

In summary, we investigated the properties of a para- 
magnetic S-spin system under Tsallis generalized statis- 
tics, on the normalized formalism. For q<0.5, a series of 
kinks appears in magnetization and entropy. This effect 
is a direct consequence of a peculiar occupation proba- 
bility, as a result of Tsallis' cut-off, where the positive 
energy states (m s < 0) vanish sharply. Additionally, the 
negative energy states (m s > 0) share non-zero occupa- 
tion probability, preventing a fully magnetized state and 
the saturation magnetization decreases with decreasing 
q. We present evidences based on experimental results 
of Amaral et al. jH| and Hebert et al.j2^] which add and 
support our previous publication J2l]], where manganites 
are suggested to be magnetically non-extensive objects. 
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VI. CAPTIONS 

Figure 1: Generalized Brillouin function for several q 
values and S=5/2, as a function of x* . 

Figure 2: Saturation magnetization -B^'Isat as a q 
function, for (a) half- integer and (b)integer S values. 

Figure 3: Probability of occupation of each energy 
level, for (a)(b)(c) S=5/2 and (d) S=2, as a function 
of x* , and various values of q. 

Figure4: (a) Unormalized (Eq. [|) and (b) normalized 
(Eq. |l(]) entropy S q , for S=5/2 and different q values, as 
a function of x* 

Figure5: Inverse of susceptibility for (a) q=1.0 and 0.8, 
for S=5/2, * function. 

Figure6: (a) Temperature dependence of the in- 
verse susceptibility for the manganite La2/3Cal/3Mn03, 
above the Curie temperature ( 267 K), at low magnetic 
field. (b)Inverse of magnetization (Eq]i"g|) as a function 
of a;* -1 , for q=0.1 
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Figure 5 - M.S. Reis et al. 
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